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MULTIPLE AND PRECESSING COLLIMATED OUTFLOWS IN THE PLANETARY NEBULA IC 4634 
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ABSTRACT 

With its remarkable double-S shape, IC 4634 is an archetype of point-symmetric planetary nebulae 
(PN). In this paper, we present a detailed study of this PN using archival HST WFPC2 and ground- 
based narrow-band images to investigate its morphology, and long-slit spectroscopic observations to 
determine its kinematics and to derive its physical conditions and excitation. The data reveal new 
structural components, including a distant string of knots distributed along an arc-like feature 40"- 
60" from the center of the nebula, a skin of enhanced [O in] /Ha ratio enveloping the inner shell 
and the double-S feature, and a triple-shell structure. The spatio-kinematical study also finds an 
equatorial component of the main nebula that is kinematically independent from the bright inner S- 
shaped arc. We have investigated in detail the bow shock-like features in IC 4634 and found that their 
morphological, kinematical and emission properties are consistent with the interaction of a collimated 
outflow with surrounding material. Indeed, the morphology and kinematics of some of these features 
can be interpreted using a 3D numerical simulation of a collimated outflow precessing at a moderate, 
time-dependent velocity. Apparently, IC 4634 has experienced several episodes of point-symmetric 
ejections oriented at different directions with the outer S-shaped feature being related to an earlier 
point-symmetric ejection and the outermost arc-like string of knots being the relic of an even much 
earlier point-symmetric ejection. There is tantalizing evidence that the action of these collimated 
outflows has also taken part in the shaping of the innermost shell and inner S-shaped arc of IC 4634. 
Subject headings: ISM: planetary nebulae: general - ISM: planetary nebulae: individual: IC 4634 - 
stars: winds, outflows 



1. INTRODUCTION 

Planetary nebulae (PNe) represent the final phases 
of the stellar evolution of low- and intermediate-mass 
stars. The shape of a PN is thought to be determined by 
the interaction of the current fast , tenuous stellar wind 
(e.g.. lPatriarchi &: Perino tto 1991) of its central star with 
the previous slow, dense wind of the Asymptotic Giant 
Branch (AGB) phase. In the so-called Interacting Stella r 
Winds (ISW) model of PN formation (jKwok et al.|[l97l . 
the critical parame ter is the density distribution of the 
AGB wind (Balick ll987f k a spherical symmetric AGB 
wind will result in a spherical PN, while a polar den- 
sity gradient in the AGB wind will result in an elliptical 
or bipolar PN, depending on the degree of the density 
gradient. 

This simplified view of PN formation has been 
challenged by the overwhelming body of observations 
of PNe showing sm all-scale features (e.g., NGC7662, 
iPerinotto et al.ll2004D. collimated bipolar outflows (e.g., 
NGC7009. iFernandez. Monteiro. & SchwarzH200l . and 
point-symmetric morphologies, including multiple point - 
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symm etric bubbles (e.g., Hen 2-47 and Ml-37, Sahai 
l2000t ) and point-symmetric collimated outflows (e.g., 
NGC6884. iMiranda. Guerrero, fc Torrefied [1999). All 
these features are difficult to interpret in the framework 
of the ISW model and have inspired new theoretical 
ideas to explain the shaping of PNe. New scenarios of 
PNe formation and shaping include magnetic collima- 
tion ( Garcfa-Scgura & Lopez 2000), the action of colli- 
mated fast winds ()Lee &; Sahall [2003L or the blowing up 
of a warped circumstellar disk by the fast stellar wind 
(iRiikhorst, Mellema. fc I ckc 200^. 

iStanghellini. Corradi. fc Schward (|1993[ ) recognized a 
new morphological class of PNe characterized by the 
point-symmetry with respect to their centers. The mor- 
phology of these PNe suggests that they are shaped by 
a mechanism involving precession, most likely the ac- 
tion of precessing or rotating collimated outflows. An 
investigation of PNe with point-symmetric morphologies 
is, thus, essential to evaluate the action of precessing 
collimated outflows in the formation of PNe and to as- 
sess the origin of the collimating mechanisms. Among 
the PNe with point-symmetric morphology, IC 4634 is 
an outstanding c ase. IC 4634 is an archetypical point- 
symmetric PN as ISchwara (J1993D found that its remark- 
able double-S-shaped morphology is reflected in its kine- 
matics: each component of the two S-shaped features ex- 
pands radially with its counterpart on the other side of 
the star expanding in opposite direction. The collimated 
outflows of IC4634 are further revealed by bow shock 
features indicating their interaction with the surrounding 
medium. Therefore, both the morphology and kinemat- 
ics of IC 4634 suggest the presence of highly symmetric 
fast collimated outflows, and the precession or rotation 
of the source that produced them. 
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Since the study of IC4634 by ISchwarzl (|1993l ). sev- 
eral studies have further investigated th e kinematical 
and morphological proper ties of IC4634 (|Haiian et al.1 
119971 iToledano et al.ll2003T) . but a detailed study is lack- 
ing. In this paper, we present a spatio-kinematical study 
of IC4634 using narrow-band archival HST WFPC2 
and ground-based images and complementary long-slit 
echelle and medium-dispersion spectroscopic observa- 
tions. The observations, described in §2, have allowed us 
to investigate the morphology, kinematics, and physical 
conditions of IC 4634 and the possible shock excitation 
of the bow shock features. A thorough description of the 
different structural components of IC 4634 derived from 
these observations is given in §3, and the action of pre- 
cessing collimated outflows is examined and modeled in 
§4. The results are discussed in §5. 

2. OBSERVATIONS 
2.1. Narrow-Band Imaging Observations 

Narrow-band HST WFPC2 images of IC 4634 in the 
H/3, Ha, [N n] A6584, and [O ill] A5007 lines were re- 
trieved from the HST archive (Proposal ID 6856, PI: 
J.T. Trauger). The total integration times were 800 s 
for the H/3 image and 1,000 s for the Ha, [N n], and 
[O III] images. These images were calibrated via the 
pipeline procedure. Cosmic rays were removed by com- 
bining different exposures obtained with the same filter, 
using standard IRAF routines. The HST WFPC2-PC 
Ha, [N n] , and [O in] images of IC 4634 are displayed in 
Figure [TJ left. The H/3 image, not shown here, has been 
used in conjunction with the Ha image to investigate the 
amount of extinction towards IC 4634 and its spatial dis- 
tribution. The mosaiced HST WFPC2 Ha and [N n] 
images have been added together to produce the deep 
Ha+[N n] image shown in Fig.\^top-right. 

The HST WFPC2-PC Ha, [N III, and [O in] images of 
IC4634 have been used to derive [N n]/Ha, [O m]/Ha, 
and [N n]/[0 ill] ratio maps. In Fig. [TJ right- center and 
-bottom, the [N n]/Ha and [O in] /Ha ratio maps are 
shown, respectively. The [N n]/[0 in] ratio map is not 
presented here because it is very similar to the [N n]/Ha 
ratio map. 

Additional narrow-band images of IC 4634 were ob- 
tained at the 2.5m Nordic Optical Telescope (NOT) of 
the Roque de los Muchachos Observatory (La Palma, 
Spain) on 2005 August 1. The images were ac- 
quired through filters that included the Ha (A c =6563 A, 
FWHM=9 A, i oxp =900 s) and Ha+[N n] (A c =6562 A, 
FWHM— 46 A, i cxp =1800 s) emission lines using the 
ALFOSC (Andalucia Faint Object Spectrograph and 
Camera) camera in imaging mode. The detector was 
an E2V 2Kx2K CCD with a pixel size of 13.5/xm, pro- 
viding a plate scale of 0'.'19 pixel -1 and a field of view 
of ~6.'5. The angular resolution during the observations, 
as derived from the FWHM of stars in the field of view, 
was 0"8-0"9. 

2.2. Echelle Observations 

High-dispersion spectroscopic observations of the Ha 
and [N n] AA6548,6584 A lines of IC4634 were obtained 
using the echelle spectrograph on the 4m Blanco Tele- 
scope of the Cerro Tololo Inter-American Observatory 
(CTIO) on 2002 June 22. The spectrograph was used in 



TABLE 1 
Log of the Medium-Dispersion Spectroscopic Observations 



Date 


Slit Position 


Position Angle 


£cxp 


Number 






[arcdeg] 


[s] 


of Spectra 


2002 June 16 


Central Star 


153 


900 


1 




Central Star 


153 


300 


i 




Central Star 


153 


(50 


3 




Central Star 


153 


20 


3 


2004 May 24 


Central Star 


153 


300 


3 




Central Star 


153 


10 


3 


2004 August 


Central Star 


90 


120 


3 




Central Star 


90 


30 


1 




Central Star 


90 


20 


1 




Central Star 


90 


10 


i 




A 


90 


300 


3 




A' 


90 


600 


3 




D 


90 


30 


1 




D' 


90 


300 


3 



the long-slit mode to obtain single-order observations of 
the Ha and [N n] AA6548,6584 A lines. The unvignetted 
slit length is ^3'. The 79 line mm -1 echelle grating 
and the long-focus red camera were used, resulting in 
a reciprocal dispersion of 3.4 A mm -1 . The data were 
recorded with the SITe 2K No. 6 CCD with a pixel size 
of 24 [im. This configuration provides a spatial scale of 
0'.'26 pixel -1 and a sampling of 3.7 kms" 1 pixel -1 along 
the dispersion direction. The slit width was set to Of! 9, 
and the resultant instrumental resolution (FWHM) was 
8 kms -1 . The angular resolution, determined by the see- 
ing, was ~l'/0. The echelle observations were made with 
the slit placed at the central star of IC 4634 and oriented 
along the position angles (PA) 60° (£ cxp =1050 s), 120° 
and 180° (£ cxp =900 s), and 150° and 166° (£ cxp =1200 s). 
High-dispersion spectroscopic observations of the Ha, 
[N n] A6584 A, and [O in] A5007 A line s were also ob- 
taine d using the IACUB spectrograph (M c Kcith et al. 
I1993D on the 2.5m NOT of the Roque de los Muchachos 
Observatory (La Palma, Spain) on 2004 June 29. The 
spectrograph was used in the long-slit mode to obtain 
single-order observations of the Ha and [N n] A6584 A 
lines in the 9 th order, and of the [O in] A5007 A line in 
the 11 th order. IACUB provides a reciprocal dispersion 
of 1.74 A mm -1 in the 9 th order and 1.43 A mm -1 in 
the 11 th order. The data were recorded with a Thom- 
son CCD with a pixel size of 19 fan. This configuration 
provides a spatial scale of 0'.'139 pixel -1 , and a sampling 
along the dispersion direction of 1.5 kms -1 pixel -1 for 
the spectrum of the Ha and [N n] A6584 A lines, and 
1.6 kms -1 pixel -1 for the spectrum of the [O in] A5007 
A line. The unvignetted slit length is 40" and the slit 
width was set to 0'.'65, resulting in an instrumental res- 
olution (FWHM) of 9.5 kms -1 for the Ha and [N n] 
A6584 A lines, and 8.5 kms -1 for the [O in] A5007 A 
line. The angular resolution, determined by the seeing, 
was ^0'.'8. The echelle observations of the Ha and [N n] 
A6584 A lines was made along a slit offset 2" West and 
46'.'4 North of the central star of IC 4634 and oriented 
along PA=0° (£ cxp =1200 s). The echelle observations of 
the [O in] A5007 A line was made along a slit placed at 
the central star of IC4634 and oriented along PA=150° 
(t oxp =1800 s). 

2.3. Medium- Dispersion Spectroscopic Observations 
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Fig. 1. — (left column) HST WFPC2 narrow-band images of the central region of IC4634 in the Ha (top), [N ll] (center), and [O III] 
(bottom) emission lines, (right column) HST WFPC2 Ha+[N ll] image of IC4634 (top) showing the outermost emission found in this 
nebula, and [N Il]/Ha (center), and [O III] /Ha (bottom) ratio maps of the central region of IC 4634. The different features identified in the 
nebula and described in the text are marked on these images. The images are displayed in a logarithmic scale, while the ratio maps are 
shown in linear scale. The grayscale upper limits of the ratio maps arc 1.0 and 4.0 for the [N n]/Ha and [O in] /Ha ratio maps, respectively. 



Medium-dispersion long-slit spectroscopic observations 
of IC 4634 were obtained using the Boiler & Chivens spec- 
trograph on the 2.1m telescope of the Observatorio As- 
tronomico Nacional de San Pedro Martir (Baja Califor- 
nia, Mexico) on 2002 June 16 and 2004 May 20 and Au- 



gust 12. Multiple observations were obtained at several 
slit positions as listed in Table 1. In all cases, the 400 
line mm -1 grating was used, the slit length was 5', and 
its width was set to 150 /zm, projecting to 2" on the sky. 
The data were recorded on a SITe3 IK CCD with pixel 
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TABLE 2 
Measured and Intrinsic Line Intensity Ratios for IC 4634 







IS 




A 




A' 




D 




D' 




Line 


F 


/ 


F 


I 


F 


I 


F 


/ 


F 


I 


H I A4340 


40±2 


44 


41±3 


45 


42±5 


46 


41±3 


45 


37±6 


41 


[O m] A4363 


5.7±0.3 


6.3 


7.5±1.6 


8.3 


6.0±0.8 


6.6 


6.0±0.4 


6.6 


5.0±0.9 


5.5 


He I A4387 


1.2±0.2 


1.3 


















He I A4471 


5.2±0.3 


5.6 


5.4±1.2 


5.8 


5.2±0.7 


5.6 


5.2±0.4 


5.6 


4.6±0.6 


5.0 


N IHA4641 + O HA4649 






















+ C IVA4658 


1.1±0.2 


1.1 


















He n A4686 


0.5±0.2 


0.5 


















[Ar iv] A4714 


1.5±0.1 


1.5 










1.0±0.2 


1.0 


1.2±0.4 


1.2 


[At iv] A4740 


1.2±0.1 


1.2 










0.5±0.2 


0.5 


0.8±0.3 


0.8 


H I A4861 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


He I A4920 


1.7±0.2 


1.7 


2.0±0.6 


2.0 


1.8±0.6 


1.8 


1.9±0.3 


1.9 


1.8±0.5 


1.8 




O m] A4959 


330±10 


322 


291±8 


284 


260± 20 


254 


290± 15 


283 


285±20 


278 




O in] A5007 


995± 30 


963 


876±18 


847 


850± 25 


822 


870± 30 


811 


835±30 


808 




N i] A5199 


0.14±0.05 


0.13 














1.7±0.5 


1.6 




Fe m] A5270 


0.15±0.06 


0.14 






0.3±0.1 


0.3 












CI in] A5517 


0.40±0.06 


0.35 


0.7±0.2 


0.6 


0.6±0.2 


0.5 


0.6±0.2 


0.5 


0.4±0.2 


0.1 




CI in] A5535 


0.58±0.05 


0.51 


0.8±0.3 


0.7 


0.6±0.2 


0.5 


0.7±0.2 


0.6 


0.6±0.2 


0.5 




N ill A5755 


0.35±0.04 


0.30 


1.5±0.6 


1.3 


2.2±0.7 


1.9 


1.2±0.3 


1.0 


1.1±0.4 


0.9 


C iv A5801-12 


0.55±0.06 


0.46 


















He I A5876 


17±1 


14.1 


16.9±1.0 


14.0 


16.2±1.2 


13.4 


17.0±1.5 


14.1 


19.7±1.3 


16.4 


[O i] A6300 


0.8±0.2 


0.6 


9.6±0.5 


7.4 


9.6±0.5 


7.5 


4.4±0.5 


3.4 


4.9±0.7 


3.8 


[S m] A6312 


1.3±0.15 


1.0 


2.3±0.5 


1.8 


2.7±0.4 


2.1 


1.9±0.4 


1.5 


2.3±0.6 


1.8 


[O i] A6363 


0.30±0.10 


0.23 


3.5±0.4 


2.7 


3.3±0.4 


2.5 


1.7±0.4 


1.3 


1.7±0.5 


1.3 


[N n] A6548 


4.2±0.4 


3.2 


37±4 


28 


42±4 


31 


20±2 


15.0 


22±5 


16.5 


H I A6563 


355±20 


264 


379±10 


283 


365±10 


272 


390±10 


290 


373±8 


278 


[N n] A6584 


11.2±0.8 


8.4 


110±6 


82 


107±6 


80 


56±5 


42 


65±6 


49 


He I A6678 


5.2±0.4 


3.8 


5.9±0.5 


4.3 


5.1±0.5 


3.8 


5.0±0.3 


3.7 


5.8±0.5 


4.3 


[S n] A6717 


0.70±0.15 


0.5 


11.2±1.5 


8.2 


10.7±1.5 


7.8 


4.6±1.0 


3.4 


4.8±1.1 


3.5 


[S ii] A6731 


1.4±0.2 


1.0 


17.5±1.5 


12.7 


14.0±1.5 


10.2 


7.8±1.2 


5.7 


8.7±1.2 


6.4 


He I A7065 


8.0±0.7 


5.6 


5.3±0.3 


3.7 


5.2±0.4 


3.6 


5.5±0.6 


3.9 


7.2±0.7 


5.1 


[Ar m] A7135 


13.5±0.8 


9.4 


18.1±1.0 


12.5 


16.5±1.1 


11.5 


18.3±1.3 


12.7 


20.2±1.5 


14.0 


He I A7281 


1.1±0.2 


0.8 










1.1±0.2 


0.8 


1.1±0.2 


0.8 


[O n] A7320 


3.0±0.3 


2.0 


9.4±0.8 


6.3 


9.0±0.7 


6.1 






9.8±0.7 


6.7 


[O II] A7330 


2.8±0.3 


1.9 


7.4±0.9 


5.0 


8.1±0.8 


5.5 






8.2±0.8 


5.6 


F(H/3) (ergs cm" 2 s" 1 ) 


1.4xl0 -11 




4.0xl0" 12 




3.5xl0 -12 




2.6xl0 -12 




2.1xl0~ 12 





size of 24/im. This configuration provides a spatial scale 
of 1//05 pixel -1 and a spectral scale of 3 A pixel -1 with 
a spectral coverage from 4240 A to 7310 A. The long-slit 
spectra were reduced and calibrated following standard 
procedures using the XVISTA package 8 . For the wave- 
length calibration, we used a He-Ar lamp. The spectral 
resolution was ~6.8 A. For the flux calibration, several 
standard stars were observed each night. 

3. STRUCTURAL COMPONENTS OF IC 4634 

The HST narrow-band images and ratio maps of 
IC 4634 shown in Fig. Q] reveal a wealth of structural 
components in this nebula. In Fig. [IJ we have labeled 
the most prominent structures in IC 4634: 

• a bright inner shell (IS), surrounded by a dim outer 
shell (OS), 

• an inner S-shaped feature formed by two knotty 
point-symmetric arcs (D-D 1 ) that present distinct 
properties along the minor axis of IS (Dm-Dm!), 

• an outer S-shaped feature formed by two pairs of 
bow shock-like structures (A- A' and B-B') along 

8 XVISTA was originally developed as Lick Observa- 
tory Vista. It is currently maintained by Jon Holtz- 
man at New M exico State University and is available at 
http://ganymede.nmsu.edu/holtz/xvista 



different PAs and located at different distances 
from the central star, 

• a pair of bow shock- like structures, C-C , along 
PA-120 , 

• a skin of enhanced [O in] /Ha ratio (S) that en- 
closes the inner shell and bow shock-like structures, 

• an outer envelope (Env) surrounding all previous 
components and where some individual features 
(F-F and G) can be distinguished, 

• a string of knots (H) distributed along an arc-like 
feature at ~40" North of IC 4634 and accompanied 
by diffuse emission (Hd). 

In the following, we will describe in greater detail the 
morphology, kinematics, excitation, and physical condi- 
tions of the different structural components present in 
IC4634. 

The H/3 fluxes, and measured and intrinsic line inten- 
sity ratios of the emission lines detected in the medium- 
dispersion spectra in regions IS, A- A', and D-D 1 are 
listed in Table 2. The Ha to H/3 ratio measured in 
the medium dispersion spectra and in the HST images 
imply an extinction coefficient, ch/3 = 0.34. The HST 
Ha/H/3 ratio map indicates that there are no signif- 
icant extinction variations across IC4634. The mea- 
sured line intensity ratios have been dereddened ac- 
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TABLE 3 
Plasma Diagnostics and Ionic and Elemental Abundances for IC 4634 



IS 



Plasma Diagnostics 

N c [S II] (cm- 3 ) 
T e [N II] (K) 
T c [O in] (K) 



10,600 

13,150 

9,850 



2,800 

9,600 

11,350 



1,500 
11,600 
10,700 



3,800 
11,550 
10,550 



5,200 
10,250 
10,050 
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He/H 
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O/H (xlO 4 ) 
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N/H (xlO 4 ) 
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Fig. 2. — HST WFPC2 narrow-band images of the inner shell of IC4634 in the Ha (left), [O III] (center), and [N II] (right) emission 
lines. 



cordin gly using the wav el ength -dependent extinction law 
from ISavage fc Mathisl (|1979f ). The plasma diagnos- 
tics (T e and N e ), and ionic and elemental abundances 
listed in Table 3 were derived from the intrinsic line in- 
tensity ratios listed in Tab. 2 for the different regions 
of IC4634 usin g the nebular abundance package ELSA 
(| Johnson et al.ll2006f) . 

3.1. The Inner Shells of IC 4634 

From the center outward, the first component of 
IC4634 is the inner shell, IS, around the central star. 
This inner shell has a size of 2'.' 5x5'.' 2 and a major axis 
along PA~150°. The edge of the inner shell is sharp 
and its surface brightness is enhanced with respect to 
the shell interior. This morphology suggests a thin shell, 
although it must be noted that an intricate pattern of fil- 
aments is superimposed on the shell. The shell brightness 
is also notably enhanced along its minor axis, marking 
a bright equatorial belt. The shape of the inner shell 
can be roughly described as elliptical (Figure [2]), but, at 



the tips of its major axis, several protrusions of differ- 
ent sizes extend outwards. Further distortion of the shell 
shape is introduced by the bending of the shell edge at 
its equator, making the edge of the inner shell along its 
major axis to show a subtle S-shape. It is worthwhile to 
note that the central star of IC 4634 is not located at the 
exact center of the shell, but displaced by ~0'.'20 towards 
the South. 

The echellogram of the [O in] emission line along the 
major axis of the inner shell at PA 150° (Figure [3]) re- 
veals a prolate expanding shell with an expansion veloc- 
ity at the position of the central star of 15.3 km s . 
The systemic radial velocity in the Local Standard of 
Rest (LSR) system is -20.1 km s -1 , in a greement with 
other estimates (e.g.. iDurand et al] [l998). The line tilt 
is small, indicating a low inclination of the shell with 
respect to the plane of the sky, but the S-shaped pat- 
tern of the line suggests accelerating features at the tips 
of the shell major axis. These blowout structures cor- 
respond to the protrusions observed at the tips of the 
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hel sys 



Fig. 3. — Echellograms of the Ho, [N II] A6584 A, and [O III] A5007 A emission lines along PA 150°, and echellogram of [N n] A6584 A 
emission line and [N Il]/Ha line ratio map along PA 60° of the inner shell of IC4634. In the [N II] echellograms, the emission is dominated 
by the adjacent regions D-D' and Dra-Dm' , but still a hollow cavity can be distinguished, as marked by an ellipse in the [N II] echellogram 
along PA 60°. The echellograms are displayed in a logarithmic scale, while the line ratio map is shown in linear scale with lower limit 
and upper limit 0.25. 



major axis of the shell, with the Northwest protrusion 
receding from the observer and the Southeast protrusion 
approaching. The echellogram of the [N n] emission line 
along the major (PA=150°) and minor (PA=60°) axes 
(Fig. [3]) also reveal an expanding shell, but its kinemat- 
ics cannot be studied in detail in this line because the 
bright emission from the adjacent D-D 1 and Dm- Dm' re- 
gions overwhelms the emission from IS. The expansion 
velocity at the position of the central star is illustrated in 
the [N ii] echellogram along the minor axis by an ellipse 
with spatial semi-axis ~l'/3 and expansion velocity of 18 
km s _1 . In the Ha echellogram along the major axis of 
IC 4634, also shown in Fig. [3J the thermal broadening of 
the line is larger than the line split and no shell can be 
clearly identified. 

The inner shell of IC4634 is the region where the 
[N Il]/Ha ratio map (Fig. [T]) shows its lowest values, 
~0.04. The high excitation is confirmed by the [N Il]/Ha 
line ratio map presented in Fig. [3] that shows this ratio 
to be nearly zero in the central regions of IC 4634. The 
[O hi] /Ha ratio (Fig. [I]) is very flat in this region, with a 
typical dereddened value of ~3.0. This is also the only re- 
gion where emission lines of He n are detected in the low- 
dispersion spectra (Tab. 1), although the intensities of 
these lines are low with /(He n 4686 A)//(H/3) ~ 0.005. 
Besides the [Ar iv] AA4715,4740 A lines, no other high- 
excitation lines are present in the spectrum of the inner 
shell of IC 4634. The [Fe in] A5270 A emission line and 
the Wolf-Rayet complex of C iv A5801-12 A are weakly 
detected in this spectrum. 

The electron density and temperature of the inner shell 
of IC 4634 can be derived using the density sensitive [S ii] 
AA6716,6731 A doublet line ratio and the temperature 
sensitive [N ii] and [O in] auroral to nebular line ra- 
tios. Using the line intensities listed in Tab. 2, we have 
computed an electron density of ~10,600 cm -3 from the 



[S ii] line ratios, and electron temperatures of 13,000 K 
from the [N ii] lines and of 10,000 K from the [O III] 
lines. These values are consistent with previous esti- 
mates of the physical conditions in the central region of 
IC46 34 (lAller fe Czvzak1ll983t lde~Freitas Pacheco et al.l 
Il992t iHvung et al.l Il999h that also find higher electron 
temperature from the [N ii] lines than from the [O in] 
lines. The nebula does not show any evidence of N or 
He enrichment, thus it can be clas sified as a type II PN 
based on its chemical abundances (|PeimbertJll978f ). 

The inner shell of IC 4634 is surrounded by a faint shell 
of radius 3" best seen at PAs 85°-150° and 265°-330° 
(Fig. [2]), the outer shell OS, as marked in Fig. [TJ At 
other PAs, the low surface brightness of this shell is over- 
whelmed by the bright D-D 1 and Dm-Dm! features. The 
visible portion of OS suggests that it is round or, at least, 
less aspherical than the innermost IS. 

3.2. The Inner Point- Symmetric Structure 

One of the most prominent features in IC 4634 is the 
pair of low-excitation arcs that we named D-D 1 in Fig.[TJ 
A close-up of these features is presented in Figure 01 D- 
D 1 have a point-symmetric morphology, bending around 
IS in arc-like features that extend from the IS equa- 
tor, at PA ~60° (240°), to the IS tips, at PA ~330° 
(150°). These features do not look like filaments, but 
rather like broad (~ 1"), clumpy strips, particularly in 
the [N n] image. The morphology is reminis cent of the 
low-excitation pola r caps of NGC 6543 (Mir anda fe SolJ 
1992; Balick 2004) and of the south fuzz of NGC 2392 
(jO'Dell et al.ll2002ft . 

The surface brightness of these structures does not 
change significantly in the [N n] line, but there is a clear 
enhancement of the Ha and [O in] emission in the equa- 
torial regions that we denote as Dm-Dm! . To investigate 
these changes in more detail, we have used the well estab- 
lished [O in] /Ha vs. [N n]/Ha line diagnostic diagram. 



Multiple Collimated Outflows in IC 4634 



!*™l 


30 


s 


* 

32 


a 




2 


34" 


2 




a 




a 


* 

36 


o 




a> 




u 


-31*49'38 







ap 


: Ha 


I [n n] g^fe 

-34 


i^l-jSj/Ha ^ 


~ 


+27 




%] 


: +21 "J 

+17 ' 

. I I 1 1 I 1 1 1 1 I 1 1 1 1 I I 1 1 M 


"i 
■ i ■ -1 i i ■ ■ 1 ■ 1 ■ m ' "-!-'■ "■■"!'-■ '..i * t- f t. 



lT^Ol^STS 33?6 33?4 33?2 33?6 33?4 33f2 

Right Ascension (J2000) 



3376 33*4 3372 



Fig. 4. — HST WFPC2 narrow-band images of IC4634 in the Ha (left) and [N n] (center) emission lines, and [N Il]/Ha ratio map 
(right) emphasizing the emission from regions D-D' and Dm-Dm! . The radial velocities with respect to the systemic velocity along the 
D-D 1 arc-like features and at the location of Dm-Dm' are marked on the [N II] image. The images and ratio maps are presented in a linear 
scale. The greyscale of the ratio map has a lower limit of and and an upper limit of 1. 
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Fig. 5. — (left) [O m]/Ha vs. [N Il]/Ha diagnostic diagrams of different morphological components of IC 4634. The exact spatial regions 
from which these diagnostic diagrams have been extracted are overlaid on an Ha image of IC4634 (right). 



To construct this diagram, we have selected the rectan- 
gular regions encompassing D-D 1 and Dm-Dm' shown in 
Figure fright and computed the dereddened [O in] /Ha 
and [N n]/Ha ratios at every pixel with surface bright- 
ness above a threshold value of 3a. The resulting val- 
ues are plotted in the [O in] /Ha vs. [N n]/Ha diagrams 
shown in Fig. [SJ-Ze/it. In these diagrams, D and D 1 present 
very similar distributions in their excitation conditions, 
and the same applies between Dm and Dm! . These show 
slightly higher [O m] /Ha values than D-D ', but very 
different [N n]/Ha distributions, with Dm-Dm' values 
lower than ~0.15, while D-D 1 present values of this ra- 
tio up to 0.6. Therefore, there is a notable segregation 
in the [O III] /Ha vs. [N n]/Ha plane between D-D 1 and 
Dm-Dm! ', indicating different excitation conditions. 



The different nature of D-D' and Dm-Dm! is further 
substantiated by their kinematics in the [N n] echel- 
logram at PA 60° (Fig[3]). In this echellogram, four 
distinct compact knots are detected at the location of 
D-D 1 and Dm-Dm! . Of these four knots, the pair of 
knots located at the same side of the nebula have op- 
posite velocities, i.e., D and Dm have very similar loca- 
tions, but very different radial velocities, and the same 
applies to D 1 and Dm! . The position- velocity (PV) map 
of the [N n]/Ha ratio (Fig. [3]) shows that the brighter 
two [N II] knots, corresponding to D-D 1 , have lower ex- 
citation, while the two fainter knots show higher excita- 
tion and can be identified with Dm-Dm!. Therefore, D 
and Dm! are blueshifted, with systemic radial velocities 
of -21 km s _1 and -27 km s _1 , respectively, while D 1 
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and Dm are redshifted, with systemic radial velocities of 
+21 km s _1 and +27 km s _1 , respectively. 

The kinematics along D-D 1 has been further mapped 
at PAs 150°, 166° and 180°. In the corresponding echel- 
lograms, D-D 1 appear as bright, compact knots with a 
velocity FWHM of -18 km s" 1 in the [N n] A6584 A 
line. D is blueshifted and its radial velocity shows little 
variations along this structure as shown in Fig. |4j D' 
is redshifted and shows velocity variations symmetric to 
those presented by D. It must be noted that these ve- 
locities are opposed to the line tilt shown by the inner 
shell, whose Northwest protrusion is redshifted, while its 
Southeast tip is blueshifted. 

The physical conditions at the tips of D-D 1 have been 
derived using temperature and density sensitive line ra- 
tios obtained from the medium dispersion spectroscopy. 
Using the line intensities listed in Tab. 2, we have com- 
puted an electron density of 3,800 cm -3 for D and 
5,200 cm" 3 for D 1 , i.e., the density of D-D 1 is 2-3 times 
lower than the density of the inner shell. The [N n] elec- 
tron temperature is found to be 11,500 K and 10,200 K 
for D and D' , respectively. The [O in] electron temper- 
ature is 10,500 K for D and 10,000 K for D 1 . 

3.3. The System of Bow Shock-Like Structures 

The system of bow shock-like structures observed at 
PAs -150° (A- A'), -166° {B-B), and ~120° (C-C) is 
one of the most remarkable features of IC 4634. There 
are very few exa mples of bow shock structures among 
PNe (e.g., IC4593. ICorradi et al.|[l996T) . In the following 
subsections, we describe in further detail the structure, 
excitation, physical conditions, and kinematics of these 
features. 

3.3.1. Small-Scale Structure 

The Ha, [N ii] and [O in] composite pictures of A 
and A' shown in Figure [5] reveal a wealth of small-scale 
structures in these features. The bow shock-like feature 
A lies 10" NW of the central star. The counter bow shock 
(feature A') is precisely symmetric on the opposite side 
of the central star (see Fig. [1]). A chain of knots extends 
towards the SE from structure A and towards the NW 
from structure A', ending in a smaller sized bow shock 
labeled B and B' in the NW and SE, respectively, that 
lie at —7" from IC 4634 central star. 

Both in A and A' there are well-defined bow-shaped 
structures with sharp edges that emit predominantly in 
[O in] (shown in blue in Fig. [BJ- top). Fainter [O in], [N n] 
and Ha diffuse emission is observed in the region between 
the bow shock and the inner nebula. In the brightest re- 
gions of A and A', several knots can be resolved appear- 
ing distinct in [N n] (Fig. O, but surrounded by diffuse 
emission in [O in]. The twisted distribution of these [N n] 
knots is indicative of the development of instabilities at 
this location. 

The bow shock structures B and B' show faint wings 
that are also more prominent and extended in the [O in] 
image. Only one knot, located at the cap of the bow 
shock structure, is resolved within B. The [O in] image 
shows that the spatial distribution of this knot is con- 
siderably more extended than the corresponding Ha and 
[N n] intensity peaks. In the emission line images of the 
counter bow shock B' , three knots are resolved. 



C and C (Fig. [1} are a pair of low-emissivity features 
located along an axis at PA —120°. In the [N n] image, 
they show pointed, triangular tips, that are surrounded 
by round, cap-like structures in the [O in] and Ha im- 
ages. At the outermost skin of these cap-like structures, 
the [O in] /Ha line ratio is enhanced. 

3.3.2. Excitation and Physical Conditions 

The [O in] /Ha and [N n] /Ha ratio maps show abrupt 
changes in these line ratios throughout A- A', B-B 1 , and 
C-C (Fig. 0. At the leading edge of the bow shock struc- 
tures, the [O III] emission is enhanced, with [O in] /Ha 
values of 2-3. At the brightest regions of A- A' and B-B' 
(i.e., at the emitting knots of the bow shock structure), 
the [O in] /Ha ratio declines to 1-1.5, while the [N n] 
emission is enhanced, with [N n]/Ha values of 0.8-0.9, 
significantly raised from the value of 0.05-0.15 behind the 
bow shock. Similarly, there are abrupt changes in these 
ratios across the B-B features (Fig. [6]), with the lowest 
[O in] /Ha ratios and the highest [N n]/Ha ratios at the 
location of the knots (i.e., at the head of the bow shock). 
The spatial offset between the [O in] enhanced cap- 
shaped bow shocks and the low-ionization [N n] bright 
knots is ~0'/3. We note that the wings at the A-A' bow 
shocks show the highest excitation among these struc- 
tures, with [O in] /Ha values of ~4, and [N n]/Ha values 
lower than —0.2. 

As in §3.2, we have computed the [O m]/Ha vs. 
[N n]/Ha diagrams for the different rectangular regions 
encompassing A-A', B-B', and C-C shown in Fig. [5} 
right. In these diagrams, shown in Fig. [Sj-Ze/t, it is no- 
ticeable the strong similarity between the distribution of 
the data points of each structure in the NW and SE 
regions. Data points from structures A-A' and B-B' 
tend to occupy an extended region in Fig. ®left, imply- 
ing large variations in the excitation conditions on small 
scales. Most points fall in the region of the diagram with 
[O in] /Ha values ~2-3 and [N n]/Ha -0.05-0.20, cor- 
responding to the gas behind the bow shock. There is 
a locus in this diagram characterized by low values of 
[O in] /Ha and large values of [N n]/Ha, which corre- 
spond to the bright emitting knots. On the other hand, 
there are also points characterized by large values of the 
[O in]/Ha ratio, >3, and very low values of the [N n]/Ha 
ratio; these points correspond to the "skin" of [O in] ob- 
served at the leading edge of the knots. The excitation 
of the bow shock-like structures B-B' is similar to that of 
A-A' . Most of the data points of B-B fall in the region 
having [O in] /Ha in the range 2-3 and [N n]/Ha ranging 
from 0.05 to 0.20. The points with low [O in] /Ha val- 
ues and high [N n]/Ha values corresponds to the [N Ii] 
bright knots. Note the presence of several data points 
with large [O in] /Ha values (ranging from 3 to 4), which 
correspond to the leading edge of the bow shocks. 

The data points from features C-C occupy a small 
region in the [O in] /Ha vs. [N n]/Ha diagram (Fig. [SJ- 
right). Most points fall in the region with large [O in] /Ha 
ratio (~ 2.25 to 3.5) and low [N n]/Ha ratio (< 0.15). 
Data points corresponding to the outer "skin" of [O in] 
show [O in]/Ha larger than 3.5. There are noticeable 
differences between the distribution of the data points of 
A-A' and B-B' in the [O III] /Ha vs. [N Il]/Ha diagram, 
and this of C-C. In particular, the data points of C-C 
have a smaller scatter in the [O in] /Ha vs. [N n]/Ha dia- 
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gram than those of A-A' and B-B 1 , implying small varia- 
tions in the excitation conditions across C-C '. Moreover, 
C-C have higher excitation than A- A 1 and B-B' . 

Our ground-based medium-dispersion spectroscopy 
(Tab. 2) cannot resolve the rapidly changing ionization 
structure seen in the HST images, but the use of addi- 
tional line diagnostic diagrams from a variety of emission 
line ratios can provide further information on the exci- 
tation and physical conditions of A- A'. To better com- 
pare the overall excitation of A- A 1 with that of the inner 
shell, IS, and D-D , the line diagnostic diagrams shown 
in Figure [7] also include emission line ratios from these 



regions. The spectra of A-A 1 show enhanced emission 
of low-excitation lines: the relative intensities of [O i], 
[O n], [N n], and [S II] are up to 10 times higher than 
those in the spectrum of IS, while the relative intensities 
of [O i] , [N II] , and [S n] are larger than these of D-D by 
a factor of ~2. The relative intensity of the [O in] lines 
of both A-A' and D-D is only 10%-15% lower than these 
in 75, while the He I and [Ar in] lines have similar inten- 
sities among these different regions, indicating that the 
radiation field reaching both A-A' and D-D , modified 
by the absorption in the IS, is similar in both cases and 
capable of ionizing species with ionization potential < 30 
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Fig. 7. — Different diagnostic diagrams for the inner shell IS 
(solid dots), the arc-like features D-D' (open squares), and the bow 
shock-like features A-A' (stars) as derived from the undereddened 
line ratios measured in the low-dispersion spectra. 



eV. We have obtained the [O m]/[0 n] intensity ratio, 
which is indicative of the ionizing parameter for a pho- 
toionized gas, for the three different regions. The largest 
value (i.e., the largest ionizing parameter) corresponds 
to the inner shell, as expected, with a value of ^250. A- 
A' and D-D 1 show lower values in the narrow range from 
~70 to ~80, thus indicating that these structures have 
similar local ionizing conditions (somewhat lower in A- 
A'). Even though the ionizing parameters are similar in 
D-D 1 and A- A', the latter shows stronger emission of low 
excitation lines (see above), indicating shock excitation. 
Fig. [7] also reveals the progressive change in the [N n] 
and [O in] temperatures with distance from the central 
star of IC 4634: the [N n] 6584/5755 emission line ra- 
tio increases (and so the temperature inferred from this 
line ratio decreases) from IS to A-A', while the [O in] 
5007/4363 emission line ratio decreases (and so the tem- 
perature derived from this line ratio increases) as we 
move outwards from regions of low [N li]/Ha values at IS 
to regions with larger values of this ratio. The electron 
temperatures of A and A' are, respectively, 10,200 K and 
12,200 K from the [N n] lines, and 10,800 K and 11,400 K 
from the [O in] lines. Similarly, the electron density from 
the [S n] doublet line ratio is found to be 2,800 cm~ 3 in 
A and 1,500 cm~ 3 in A', i.e., 3-7 times lower than in IS. 

3.3.3. Kinematics 

The large opening angle of these bow shock-like struc- 
tures (A and A' have almost flat morphologies) suggests 
a low inclination angle with the plane of the sky. This 
is confirmed by the PV diagrams of the bow shock-like 
structures A-A' and B-B' displayed in Figure[8]that show 
low radial velocities with respect to the systemic veloc- 
ity: ±20 km s- 1 for A-A', and ±3 km s" 1 for B-B'. The 
Ha and [N n] A6584 A echellogram along PA=120° (not 



shown in Fig. [8]) shows faint emission from C and C 
moving at radial velocities with respect to the systemic 
velocity of ±30 km s _1 . We note that A, B, and C (A', 
B' , and C) are receding from (approaching) us, opposite 
to the motion of D (D 1 ). 

It is also interesting to note the detailed line shape in 
the PV diagrams of Fig. M at A-A' and B-B'. At the 
location of A-A' , the Ha and [O in] emission lines are 
wedge-shaped with the widest side at the line tip (the 
leading edge of the bow shock), while the [N n] line shape 
is nearly round. The FWHM of the Ha emission line of 
both A-A' increases from ~30 km s _1 at the narrow end 
of this wedge-shaped feature to ~70 km s _1 at its wide 
end. A close examination of the Ha and [O in] 5007 A PV 
diagrams reveals a velocity gradient at the leading edge of 
A-A', whereas the velocities in the region behind the bow 
shock seem to be rather constant. The lowest velocities 
are found at the leading edge of the bow shock, facing 
away from the central star of IC 4634, where the radial 
velocity of the Ha line decreases ~10 km s _1 within 1". 
Across B-B' , the radial velocity is roughly constant. The 
Ha emission line profiles show a mean FWHM of ^40 km 
s _1 . The emission ends in a sharp drop at the outer edge 
of B-B' . 

Finally, the string of [N n] knots that arises from the 
edge of the A and A' bow shocks may be suggestive of 
a trail of material left in their motion forward. Indeed, 
the radial velocity in the outer low-ionization arcs shows 
an abrupt change from the bow shock structure to the 
linear string of knots, suggesting a deceleration. We will 
show in §4.2 that the interaction between a precessing 
collimated outflow with time-dependent velocity and the 
surrounding medium offers an alternative explanation to 
the string of [N n] knots. 

3.4. The [O in/ Skin 

It has been mentioned that the [O in] /Ha ratio is 
clearly enhanced in cap-like structures just outside the 
bow shocks A-A', B-B', and C-O . The [O in] /Ha ra- 
tio map (Fig. \\\ right- center) reveals that these [O in] 
enhanced cap-like structures extend on a strip that in- 
scribe the inner regions of IC 4634, like a thin shell or 
skin of enhanced [O in] emission. Only the outer enve- 
lope (Env) and helical structures (H and Hd) are located 
outside this skin of enhanced [O in] /Ha ratio. 

The presence of a skin of enhanced [O in] /Ha ratio 
is rare among PNe; only in NGC 6543 a structure of 
this kind has been reported (|Balickl r2004f ). In his de- 
tailed analysis of narrow-band HST WFPC2 images of 
NGC 6543, Balick noticed that there is no observational 
artifact (contamination of [N n] emission in the Ha im- 
age or incorrect correction for wavelength-dependent ge- 
ometric distortions in the camera optics), nor ionization 
effects or local variations of the O/H abundances or the 
electronic temperature in NGC 6543 able to produce the 
observed increase of the [O in] /Ha ratio at the nebular 
edge. Observational artifacts can also be dismissed as 
the origin of the observed morphology in the [O in] /Ha 
ratio map of IC 4634, as the emission in [N n] is much 
weaker than in Ha, as well as a local variation of the 
O/H abundances which are rather constant throughout 
the nebula. Emission from the main nebula scattered in 
the PN outer layers can also be excluded as the [O in] /Ha 
ratio in this skin is significantly different from the ratio 
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Fig. 8.— Echellograms of the Ha, [N n] A6584 A, and [O in] A5007 A emission lines along PA 150° 
[N n] A6584 A emission lines along PA 166° of the bow shock-like features A- A' and B-B' . 



and echellogram of the Ha and 



value measured in the bright, innermost regions of the 
nebula. The only plausible cause of the observed en- 
hancement in the [O in] /Ha ratio is a local increase of 
0++/H+ or T c . 

The origin of the skin of enhanced [O h i] /Ha ratio i n 
NGC6543 was concluded to be uncertain (|Bahckll2004h . 
In IC4634, this structure is related to the bow shock 
structures A-A', B-B , and C-C and, therefore, it can be 
speculated that it originates in the interaction of fast col- 
limated outflows with the nebular material. The shocks 
produced by fast collimated outflows would excavate a 
cavity in the low density nebular envelope and propagate 
outwards, inducing a marginal increase of the electronic 
temperature in a forward shock that enhances the [O in] 
as also observed in w ind-blown bubbles around WR stars 
(jGruendl et al.ll2000| ). A similar origin can be attributed 
to such structure in NGC 6543, as suggested by the caps 
of enhanced [O in] /Ha associated with the jet- like fea- 
tures of NGC 6543 (see Figure 2 of Balick 2004). 

3.5. The Outer Envelope 

The HST images of IC 4634 reveal an envelope of faint 
emission surrounding its central regions that we have 
labeled as Env in Fig. [TJ This envelope has a patchy 
appearance, with a distinct arc towards the East, and 
individual filamentary (F-F) and arc-like (G) features 
towards the West. Some of these features are best seen 
in the [O in] /Ha ratio map shown in ¥ig.\\^right-bottom. 
Overall, the morphology of the envelope can be classified 
as elliptical, with its major axis along PA~120°. 

Since the envelope is more clearly detected in the Ha 
and [O ill] images than in the [N n] image, it suggests 
a high excitation. Its surface brightness is low, up to 
~1,000 lower than that of the inner shell in the Ha line. 
If we assume the same electron temperature as in the 
inner shell, its density can be scaled down from that of 



the inner shell to wlOO cm -3 . 

Finally, the envelope is detected in the high-dispersion 
spectra, especially in the echellograms of the Ha line. 
Its emission shows a broad, unresolved line with a radial 
velocity similar to the systemic velocity. This structure 
seems to be an inert irregularly shaped, low-density outer 
envelope. 

3.6. The Outer Helical Structure 

The existence of a distant string of fai nt knots associ- 
ated w ith IC 4634 was first reported by iGuerrero et al.l 
(2004). This feature, labeled H and Ed in Fig.^right- 
top, is shown in greater detail in the deep NOT Ha+[N n] 
image presented in Figure [9) In this figure, H and Hd 
are composed of several knots detected from PA ~— 5° 
to 20° at angular distances 37"-64" from the central star 
of IC 4634. In the HST images, the knots in the struc- 
ture H are resolved into compact cores surrounded by 
clumpy emission or show the appearance of a string of 
knots. Morphologically, component H forms an arc-like 
string of knots whose orientation is different from that of 
the major nebular axis, but similar to this of components 
B-B 1 and D-Lf . Component Hd consists of several faint 
diffuse knots, including a large (~ 25") corkscrew-shaped 
structure oriented along the South-North direction. We 
have obtained deep narrow-band images to search for 
a possible Southern counterpart of H and Hd, but our 
search has yielded negative results. 

The Ha and [N n] emission from H and Hd is detected 
in the high-dispersion spectroscopic observations. The 
emission lines can be fit with a single Gaussian profile 
with FWHM of -25 km s^ 1 in the Ha line and -20 
km s _1 in the [N n] A6584 line. The systemic radial 
velocity derived from these lines is very similar to the 
nebula systemic radial velocity. 
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Fig. 9. — Grey scale representation of the NOT Hct+[N III image 
of the string of faint knots forming component H (see Fig.[lJ. The 
arrows and closed curves overlaid on the image mark the nebulosi- 
ties associated with this component, while the white stars indicate 
the position of field stars overimposed onto some of the line emis- 
sion knots, as identified by comparing this narrow-band image with 
an HST WFPC2 broad-band V image. The origin (0,0) is at the 
position of the central star of IC 4634. 



4. A PRECESSING COLLIMATED OUTFLOW IN 
IC 4634 

4.1. Observational Evidence and Similarities with 
Herbig-Haro Objects 

The morphology of the bow-shaped features A-A', 
B-B' , and C-C resembles that of Hcrbig-Haro (HH) 
jets. These jets are morphologically characterized by 
chains of aligned knots with bow shock-like appearance, 
and a leading bow shock, the "head" of the jet, where 
the supersonic flow slams into the surrounding mate- 
rial. Observations of the knots of HH jets have been 
interpreted by several authors as the result of time- 
depe ndent variations in the velocity of the flow (see, 
e.g., iRaga et aLlll990fh Such a jet model with variable 
ejection velocity has been shown to reproduce the over- 
all morphology and kinematical properties of the high- 
yelocity jets of the protq-plane tary nebula Hen 3-1475 
(Vela zquez. Riera fc Rag a 2004). The bow shock struc- 
ture of A-A', B-B', and C-C is consistent with those 
features being the result of a jet, ejected with variable 
velocity, interacting with the AGB remnant. 

The suggestion drawn from the morphology of IC 4634 
is reinforced by the kinematics of A-A'. Its Ha and 
[O in] PV diagrams clearly show a wedge-shaped fea- 
ture analogous to those observed in several HH objects 
(e.g- jBohm fc S olf 1985ft . A wedge-shaped feature in the 
PV diagram is expected for a bow shock moving nearly 
on the plane of the sky, as predicted by the traditional 
"3/2" -dimension bow shock model which is a reason- 
able i apjDroxirmrt^ro^ HH 
jet (|Raga fc Bohmlll985Lfl98llHartigan et al.lll987J ). In 
such case, the radial velocity dispersion increases sud- 
denly at the stagnation point of the bow shock, produc- 
ing a line broadening that looks like an horizontal edge 
in the PV diagram. The amount of broadening of the 
emission line profiles can be related to the velocity of 



the bow shock (|Hartigan et alJ ll987). In the case of A- 
A' , the observed width of the line would indicate a bow 
shock velocity of ~ 100 km s" 1 . We conclude that the 
PV diagrams of A-A' are, at least qualitatively, compat- 
ible with the predictions of a leading bow shock (i.e., the 
"head" of a jet) nearly on the plane of the sky. 

Further support for this scenario is provided by the 
spectral properties of A-A'. The dereddened emission 
line ratios in these regions with intense emission in low- 
ionization lines are reminiscent of shock-excited nebula 
such as HH objects. While shocks might be collision- 
ally exciting the low ionization emission lines of A-A' , 
a detailed comparison of the observed emission line ra- 
tios with spectra p redicted by photoionized shock models 
(e.g., Dopita 1997) is not straightforward, since the bow 
shocks in IC 4634 are illuminated by the stellar ionizing 
flux from the post-shock direction. A more realistic com- 
parison is enabled by the axisymmetric numerical simu- 
lations of a shocked, dense cloudlet moving away from a 
source of ionizing photons through a uniform and pho- 
toi onized environment th at h ave been recently o btained 
bv lRierafc Raga I (|2007ft and IRaga et"al] (|2007ft . These 
simulations not only produce synthetic spectra that are 
qualitatively consistent with the spectra of A-A', but 
also reproduce the ionization stratification of A-A' , with 
the gas in the downstream region being more highly ion- 
ized than the outward-facing edge and the [O ill] emis- 
sion showing a larger extension towards the photoioniz- 
ing source. 

It is thus tempting to interprept the point-symmetry 
morphology of IC 4634 as the result of moderate veloc- 
ity jets with a time-dependent ejection velocity. Further- 
more, the point-symmetric morphology of IC 4634, which 
gives its S-shaped appearance to the nebula, and the 
varying sign of the radial velocity are commonly inter- 
preted as the direct result of a time-dependent direction 
of ejection of the source. In the next section, we explore 
into further detail this supposition that can, in general, 
be of i nterest to other poin t- symmetric PNe (ICliffe et all 
19951: iLee fc S ahai 2 0"ol iGarcia-A rredondo fc Frank 
20041: Velazquez, Riera fc R aga 2004; Sok er fc Biskerl 
2006). 



4.2. Modeling a Precessing Collimated Outflow in 
IC4634 

4.2.1. Initial Conditions and Numerical Methods 

The morphology and kinematics derived for IC 4634 
place important constrains and clues on some of the pa- 
rameters describing the motion of such a precessing out- 
flow. The bow shock morphology of A-A', and the spatial 
distribution of the [O in] and Ha emission, preceding the 
[N 11] emission, reveals shock excitation. The observed 
radial velocities are low, suggesting that the angle be- 
tween the precession axis and the plane of the sky, the 
inclination angle i, is small. We will thus assume that 
the precession axis is on the plane of the sky. Finally, the 
semi-angle of the aperture of the precession cone, a, can 
be derived by measuring the angle subtended by A and B 
with respect to the central star. This angle is estimated 
to be ~11°. 

Numerical simulations were car- 

ried out with the 3D code Yg uazu-a 
l|Raga. Navarro-Gonzalez. fc Villagran-Munizl 2000: 



Multiple Collimated Outflows in IC 4634 



13 



iRaga et ah! |2002( ), using a 5-level binary adaptive 
grid. The dimensions of the computational domain are 
2.2xl0 17 cm, in the x- and y-directions, and4.4xl0 17 cm 
in the z-direction, i.e., 6'.'3xl2'.'7 at a distance of 2.3 
kpc, with a maximum resolution of 8.6xl0 14 cm, 
i.e., ^0'.'025. This code integrates the gas-dynamic 
eq uations by using the "flux vector splitting" scheme 
of IVan Leen (|l982f ). Together with the gas-dynamic 
equations, several rate equations for the atomic/ionic 
species are also integrated. These species are: H I, 
H ii, He i, He II, He III, C II, C in, C iv, N i, N II, 
N in, O i, O ii, O in, O iv, S ii, and S in (see details 
abou t the reaction and the cooling rates in IRaga et al.l 
|2002[) . These rate equations enable the computation of a 
non-equilibrium function for the radiative losses. Given 
a set of initial conditions for the jet and surrounding 
circumstellar medium (CSM), Yguazii-a determines 
the temperature and density distributions at a given 
time. The temperature and density distributions allow 
us to compute the emission line coefficients of the Ha, 
[N ii] A6584, and [O ill] A5007 emission lines. The 
intensity of the Ha line is computed considering the 
contributions from the recombination cascade and from 
n = 1 — * 3 collisional excitations. The intensity of 
the forbidden lines [N n] A6584 and [O in] A5007 are 
calculated by solving 5 - level atom problems, using the 
parameters of|Mendoza (|1983| ). These intensities can be 
integrated along the line of sight to produce synthetic 
emission maps and PV diagrams for a wide slit covering 
completely the working surfaces. 

A jet of number density rij is injected in the base of the 
computational domain (z = 0) at position (1.1 xlO 17 cm, 
1.1 xlO 17 cm) in the xy-plane, with a radius and length 
of 2.6xl0 15 cm and 4.3 xlO 15 cm, respectively. The jet 
density is assumed to be 10 4 cm~ 3 . The jet velocity is 
modeled by the relationship: 



vq + Av sin ( 27T — 



(1) 



where Vq is the mean velocity, Av is the amplitude of 
the velocity variation, t is the time, and r v is the period 
of the velocity variation. The symmetry axis of this jet 
precesses with period r p . 

The surrounding CSM has been assumed to be pro- 
duced by an AGB wind of constant mass-loss rate, M, 
and expansion velocity, wagb- The dependence of the 
density of such medium with the radial distance to the 
central star, r, is given by: 



p{r) 



, r ° 



(2) 



where the density, po, at radius tq is determined by the 
mass loss-rate and expansion velocity of the AGB wind as 
po = Af/(47rrg vagb)- In our simulations, the mass-loss 
rate and expansion velocity of the AGB wind were as- 
sumed to be 10 -6 Mq yr -1 and 20 km s — 1 , respectively. 
The temperature of the CSM has been set to 100 K. 

4.2.2. Results of the Numerical Simulations 

Several numerical simulations were carried out varying 
Vj, Av, r p , and r v , to produce synthetic emission maps 
and PV diagrams that can be directly compared to the 



observed images and long-slit echellograms of A- A', B- 
B' , and D-D 1 . There is a qualitative agreement between 
the synthetic and observed images and PV diagrams for 
i> =300 km s" 1 , Aw=25 km s _1 , t p =460 yr, and t v =110 
yr (Figures ITOl and [TT|) . The jet expansion velocity im- 
plies a jet mass-loss rate of ~2xlO _7 M0 yr -1 . 

Fig. [TU] displays the Ha, [N n], and [O ill] simulated 
emission maps obtained for a calculation time of 500 
yr. At this time, the jet structure achieves a length of 
3.3x10 cm, which corresponds to an angular size of 
9'.'8 at a distance to IC4634 of 2.3 kpc. The upper pan- 
els correspond to the yz-projection (i.e., the line of sight 
is along the a;- axis), while the bottom panels represent 
the a;z-projection (i.e., the line of sight is along the y- 
axis) , basically meaning the initial direction along which 
the jet is emitted. Several working surfaces are observed 
in these maps at distances of 1.4xl0 17 cm, 2.4xl0 17 cm, 
and 3.3xl0 17 cm, which are produced by the variabil- 
ity of the jet velocity, when slow gas is swept up by fast 
material. The observed morphology and spatial distri- 
bution of A- A' and B-B' are more closely reproduced by 
the yz-projection than by the ccz-projection. We notice 
that the spatial separation seen at the head of the A- A 1 
features between the Ha and [O ill] emissions are not 
reproduced by the synthetic emission maps because the 
spatial structure of the region behind the shock, where 
the [O in] emission arises mainly, is poorly resolved in 
our simulations. 

Fig. [TT1 displays the Ha, [N n], and [O in] simulated 
PV diagrams obtained at the same calculation time as 
Fig. [TUJ The centroids of the working surfaces in these 
PV diagrams do not show exactly the alternate positive 
and negative velocities observed in A- A 1 and B-B' , al- 
though they display similar sinusoidal variations with 
respect to the systemic velocity. This indicates that the 
initial direction along which the jet is emitted is not co- 
incident with neither the yz- and xz-projections, or that 
the inclination angle differs from the one that has been 
assumed. On the other hand, a close examination of the 
PV diagrams in Fig. QT] reveals a good agreement with 
the detailed kinematics of the working surfaces: the ve- 
locity widths or FWHM of the synthetic emission lines 
are comparable to the observed ones, with larger values 
of the FWHM for the Ha line, and the wedge-shaped 
tips of the Ha and [O III] emission line profiles of A- A 1 
are also well reproduced. 

5. DISCUSSION 
5.1. Multiple Shaping Agents in IC 4634 

The data presented in the previous sections reaffirm 
the known complexity of IC 4634. The nebula is com- 
posed of a series of morphological components that 
clearly reveal the action of different shaping agents. 
Among the prevalent shaping mechanisms in IC 4634, we 
shall consider the interaction of the current fast stellar 
wind with the slow AGB wind, the ionizing flux of pho- 
tons from its central star, and a series of fast collimated 
outflows. We discuss below the effects that the fast stel- 
lar wind and fast collimated outflows have had in the 
shaping of IC 4634. 

IUE observations of the central star of IC4634 un- 
covered prominent P Cygni profiles in the C in, C iv, 
N v, and O v lines superposed on the stellar contin- 
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Fig. 10. — Hq (left), [N II] (center), and [O III] (right) simulated 
emission maps of a velocity time-dependent precessing collimatcd 
outflow at a time of 500 yr. The uppers panels show the simu- 
lated emission maps for the j/ ,2 -projection, while the xz-projection 
is displayed in the bottom panels. The spatial axes show the dis- 
tance to the central star of IC 4634 in units of 10 cm. The maps 
are displayed in logarithmic scale with the flux level in units of 
ergs s _1 cm -2 sti 
each upper panel. 



vmm. These P Cygni profiles manifest the presence of 
a fast stellar wind whose terminal veloci ty has be en de- 
rived to be 3,500-4,000 km s" 1 (|Hvunget al.lfT9 99). This 
fast stellar wind has excavated a central cavity within 
the nebula, compressing the nebular material into a thin 
shell. This thin shell and the central cavity can be iden- 
tified with the inner shell IS, while the nebular envelope 
can be associated with the faint shell OS surrounding 
this inner shell. In this respect, the double-shell mor- 
phology of the inner regions of IC4634 compares well 
with other PNe (e.g., NGC6826). 

The bow shock morphology and spatial distribution 
of A- A', B-B', and C-C are highly indicative of the 
interaction of precessing collimated outflows with sur- 
rounding material. Indeed, the detailed morphology and 
kinematics of A- A' and B-B' are successfully reproduced 
assuming the interaction of a precessing jet with a time- 
dependent velocity in the range of 300 km s" 1 and a 
precession period ~460 yrs. This jet is interacting with 
material in the outer regions of IC 4634 that forms the 
envelope, Env, a region of rough elliptical symmetry that 
may represent an episode of major mass loss prior to the 
one that formed IC 4634 inner regions. In their interac- 
tion, the collimated outflows generate bow-shocks struc- 
tures and drive a forward shock that produces a high 
excitation skin. 




Fig. 11. — Simulated PV diagrams for the Ha (left), [N II] (cen- 
ter), and [O III] (right) emission lines of a velocity time-dependent 
precessing collimated outflow at a time of 500 yr. The bottom 
panels display the rrz-projection, while the upper panels show the 
yz-projection. In the panels, the horizontal axis displays the sys- 
temic radial velocity in units of kms -1 , and the vertical axis the 
distance to the central star of IC4634 in units of 10 17 cm. The 
maps are displayed in logarithmic scale with the flux level in units 
of ergs s -1 cm -2 str -1 (cm s — x ) — 1 shown in the colorbar located 
at the top of each upper panel. 

One of the new nebular components revealed in this 
work is the outermost features H and Hd. The location 
of H and Hd with respect to IC4634 and their similar 
radial velocities make it unlikely that these components 
are unrelated to the nebula. The morphology of the H 
feature, composed of a series of knots distributed along 
an arc-like structure, can be interpreted as a loop on the 
surface of an imaginary cone whose vertex is coincident 
with the central star of IC 4634. Therefore, H and Hd 
may correspond to an ancient precessing ejection. It is 
worth emphasizing that the precession axis of this an- 
cient ejection would be almost orthogonal to the axis of 
the precessing ejection that has been considered for A- A' 
and B-B' . 

There are, thus, clear signs that IC4634 innermost 
shell, IS, and the envelope, Env, formed by the action 
of the current fast stellar wind on material associated to 
episodes of mildly asymmetric mass loss produced dur- 
ing the late AGB, while the bow shock structures A- 
A' , B-B' , and C-C, and the arc-like features H and Hd 
have been produced by the action of fast, precessing colli- 
mated outflows. These collimated outflows may have also 
played some role in the shaping of the innermost regions 
of IC 4634, although the evidence is not overwhelming. 
The inner shell of IC4634, IS, is elongated and shows 
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morphological and kinematical signs of blowout along its 
tips. While the origin of the asymmetry of IS can be at- 
tributed to a density enhancement in the equator of the 
nebular envelope as suggested by the bright emission in 
the surrounding material along the minor axis in the Ha 
and [O in] images, it is worthwhile noting that the orien- 
tation and inclination of the collimated outflows associ- 
ated with A- A' are coincident with those of the blowout 
along the tips of the inner shell. Therefore, it is possi- 
ble that the collimated outflows responsible of A- A' may 
have been involved in the emergence of the blowout struc- 
tures of IS and in the development of the axisymmetry 
of this shell. This may also be the case of the inner shells 
of several PNe with collimated outflows (e.g., NGC6210 
and NGC 6884) that also show the mild point-symmetry 
exhibited by the inner shell of IC 4634. 

Similarly, the origin of the D-D 1 features may be as- 
sociated with the action of collimated outflows. The 
kinematics and morphology of D-D 1 are very similar 
to those o f the low-ionization polar caps de scribed in 
NGC 6543 (Miranda & Solf 1992; Balick 2004), although 
these structures may have different inclinations. A thor- 
ough discussion of the possible origins of k nots and low- 
ioniza tion structures in PNe is provided bv lO'Dell et al.l 
(J2002T ). As in NGC 6543, the D-D 1 arcs of IC 4634 can be 
interpreted as intrusions of low-ionization, high-density 
knots that are expanding with the s hell that surrounds 
its inner shell (jBalick fe Ha iian 2004). These features 
may represent the relic of the interaction of a collimated 
outflow with the nebular envelope. On the other hand, 
Dm- Dm', that has higher excitation than D-D 1 , and that 
are separated from these in the velocity space, can be re- 
lated to the density enhancement around the equatorial 
region of IS. 

5.2. Bow Shocks in Planetary Nebulae 

PNe display a large variety of low-ionization structures 
of which fast collimated outflows deserve special atten- 
tion because of their outstanding kinematica l properties 
(e.g.. iGoncalves. Corradi. fc Mampasoll2001[ ). The ori- 
gin and nature of the fast collimated outflows seen in 
PNe has been disputed, but the morphological, kinemat- 
ical, and physical properties of some of them seem to be 
suggestive of high-density bullet s moving supers onically 
through the nebular material ( Balick et al.lll994J ) . If this 
were the case, the interaction of the PN fast collimated 
outflows with the nebular material should lead to the for- 
mation of bow shocks whose morphologies would share 
some similarities with these of HH objects. Bow shocks 
would then be expected to be common among PNe with 
fast collimated outflows; the number of bow shocks iden- 
tified in PNe, however, is small. In IC 4634, we find one 
of the rare cases of bow shocks in PNe. The detailed mor- 
phology and kinematics of the bow shock structures A- A' 
and B-B 1 of IC 4634 can be unambiguously ascribed to a 
precessing fast collimated outflow, and the morphology 
and ionization of the bow shock structures C- C are rem- 
iniscent of the interaction of a fast, dense bullet with the 
nebular material. In view of these findings, it is worth- 
while to revisit other cases of bow shocks in PNe. 

So far, the best studied case of a bo w shock in a PN 
is that of IC4593 ([Corradi et alJll996f >. In this nebula, 
we find several systems of knots and inward facing tails 
that have been interpreted as multiple collimated out- 



flows propagating along different directions. These knots 
are preceded by outward facing caps with bow shock mor- 
phologies. The bow shock structures in IC 4593 share 
many similarities with those of IC 4634. Both have low 
radial velocities, show a decrease in the observed radial 
velocity and the broadening of the emission profile, and 
display caps of enhanced [O ill] emission as expected in 
bow shock models. All these properties suggest the mo- 
tion of medium velocity jets near the plane of the sky 
that are interacting with the surrounding material. 

The ansae of NGC 7009 are another promising example 
of fast, dense bullets ramming through the outer shells 
of the nebula. This interpretation was questioned as the 
ionization s tate of the gas de clines at the head of the 
bow shock ([Balick et al.lll998p. while the op posite is ex- 
pected in bow shocks (JHartiga n et al.lll987h . However, 
the observed ionization gradient can be reproduced by 
the numerical simulations of a shocked cloudlet moving 
away from the central star if the stellar ionizing pho- 
ton flux, which modifies the ioniz ation and excitation 
structure of the s hock, is included ( Riera fe Raga 2007; 
Raga et al]|2007D. In their study of IC 4593 bow shocks, 
Corradi et al.l ( 1996h noted that the ionizing flux from 
the central star can play a predominant role in the ion- 
ization stratification of bow shocks in PNe. At any rate, 
an examination of a deep HST WFPC2 [O in] /Ha ratio 
map of NGC 7009 reveals that the ansae of NGC 7009 
are precede d by caps of enhanc e d [O hi] as expected in 
bow shocks ([Medina et al.ll2007HGuerrero et alJl2008T > . 

Besides IC4593, IC4634, and NGC 7009, there are 
very few other cases of bow shock associat ed with PNe . 
The [O in] /Ha ratio image of NGC 6543 ([Balickl l2004h 
shows a bow shock structure enveloping the precessing- 
like coll imated outflows along the polar directions of this 
nebula (|Miranda fc SolJll992f) . although this structure is 
not discussed in detail. A bow shock structure has also 
been described in NGC 6572, although the limited spa- 
tial resolution of ground-based images do not allow the 
author s to study the ioniza tion stratification at the bow 
shock ([Miranda et al.lfl999h . A comprehensive investiga- 
tion of the occurrence of bow shocks associated with the 
expansion o f fast collimated outfl ows or bullets in PNe is 
in progress ([Guerrero et alJ l2008). 

6. SUMMARY 

Our spatio-kinematical study of IC 4634 has revealed 
new structural components in this nebula, including 
what seems to be the relic of ancient precessing colli- 
mated ejections, a triple shell morphology, and a thin 
skin of enhanced [O in] /Ha enveloping the nebula and 
the bow shock structures at the tip of the collimated out- 
flows. Furthermore, the HST images and ground-based 
echelle spectroscopy provide a detailed view of the phys- 
ical structure of the bow shock associated with a colli- 
mated outflow in a PN. Their morphology and kinematics 
have been successfully modeled using hydro-dynamical 
simulations in which a precessing fast collimated outflow 
interacts with nebular material. 

IC4634 seems to have experienced a series of mildly 
asymmetric mass loss episodes that have removed the 
stellar envelope, yielding to the current fast stellar wind. 
These periods of mass loss are interspersed with two 
episodic ejections of collimated outflows that have inter- 
acted with the nebular material, having important con- 
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sequences in the shaping of IC 4634. The most recent 
episode of ejection of collimatcd outflows is directly re- 
sponsible of the formation of the bow shock structures A- 
A', B-B , and C-C within the envelope, and has played a 
significant role in the asymmetry and orientation of the 
inner shell. The probable interaction of this collimated 
outflow with circumstellar material might be linked to 
the formation of D-D 1 . Finally, the oldest episode of 
ejection of collimated outflows formed H and Hd in the 
outermost regions. The two ejections of collimated out- 
flows took place at different times during the PN forma- 
tion. The large misalignment between the ejection axes 
of the collimatcd outflows giving raise to A- A' , B-B 1 , and 
C-C , on one hand, and H and Hd, on the other, imply 
that the collimating source had experienced important 
changes between the ejections. 

The data presented here have been taken using AL- 
FOSC, which is owned by the Instituto de Astroffsica 
de Andalucia (IAA) and operated at the Nordic Optical 
Telescope under agreement between IAA and the NBI- 
fAFG of the Astronomical Observatory of Copenhagen. 
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